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We study a pseudogap region of the mixed boson fermion system using a recent formulation 
of the renormalization group technique through the set of infinitesimal unitary transformations. 
Renormalization of fermion energies gives rise to a depletion of the low energy states (pseudogap) 
D ' for temperatures T* > T > T c which foreshadows appearance of the pairwise correlations with 

their long range phase coherence being missed. With a help of the flow equations for boson and 
fermion operators we analyze spectral weights and finite life times of these quasiparticles caused by 
interactions. 
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Many experimental measurements of e.g. the angle resolved photoemission (ARPES), tunneling spectroscopy (STM), 
specific heat as well as the optical and magnetic characteristics (see for an overview Q) point at a presence of pseudogap 
in the normal phase spectrum of underdoped cuprate superconductors. It is believed that a proper description of this 
phenomenon would be crucial for understanding of the HTSC mechanism. 

Model of itinerant fermions coupled via the charge exchange interaction to boson particles 0, for example of 
bipolaronic origin, has been shown to posses such a pseudogap in the normal phase fermion spectrum. Its occurrence 
has been obtained on a basis of the selfconsistent perturbation theory ||, diagrammatic expansion Q and also 
by solving the dynamical mean field equations for this model within the non-crossing approximation (NCA) [^). 
O \ Unfortunately, these techniques have not been able to study a superconducting phase of the model, so evolution of 
the pseudogap into the true superconducting gap could not be investigated. Only very recently we have reported 
y—{ | some prior results for the boson fermion (BF) model using the flow equation method Q where superconducting gap 
. and normal phase pseudogap could be treated on equal footing. Both gaps originate from the same source (from 
the exchange interaction between fermions and bosons) but they are differently scaled against a magnitude of this 
interaction and the total concentration of charge carriers. Still the most intriguing question concerning the evolution 
from the pseudo- to the superconducting gap is poorly understood due to technical difficulties (it is a very tough 
problem to solve the flow equations for 2 + e dimensional system) . Some attempts in this direction are currently under 
■ consideration. 

In this paper we address another important issue concerning life time effects of the quasiparticles. The high 
resolution ARPES data Q reveal that normal phase of underdoped cuprates has the marginal Fermi liquid type 
properties where the quasiparticle weight and life time are very sensitive to temperature. We want to check how the 
life times of fermions and bosons depend on temperature T and momentum using the BF model scenario. 
Flow equation for any arbitrary operator, say A(l), is given by 
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where fj is a generating operator of the continuous canonical transformation. The initial BF model Hamiltonian 
[ consists of the free part 

. S ; Ho = - Mjka + J2( E i - 2 M)&& (2) 

k,a q 

and the interaction 

H int = -t= {vk, P bl +p c ki Crt + vlph+pclfcl^ . (3) 
In order to decouple fermion from boson subsystems we have previously chosen ij in a form 
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m = - E («*.9(0«M(06 t (02*.i(oe J ,T(o - ( 4 ) 

with a kq = ej. + e q — E k+q . Using the general flow equation ([!]) we transformed the Hamilton operator H(l) so that 
in the limit I — > oo the charge exchange interaction got eliminated v k , q {oo) = 0. 

In a course of the continuous canonical transformation the initial boson and fermion operators change according to 
the following flow equations 

db(1 }} = ^<Xk,q-k(l)Vk,q-k(l)Ck,l(l)Cq-k,l(!) (5) 



dl 



E a ^wwowo< x (o- (6) 



Similar equations can be derived for their hermitean conjugates b q (l), c\ ^(l). By inspecting (||) we notice that initial 
boson annihilation operator is coupled to the electron pair operator c k ^(l)c q - k ^(l). Thereof one should in a next step 
determine a corresponding flow equation for the finite momentum pair operator, however this process will never end 
(like equations of motion for the Greens functions) . There would be the higher and higher order product of operators 
getting involved into the differential flow equations. 

Following the other works using such continuous canonical transformation ||- flicfl we terminate the flow equations 
by postulating some (physical) Ansatz. In a context of the BF model it seems natural to postulate the following 
decompositions 

b q (l) = X q {l)b q + J2 Y <i,kQ)ck,lc q -k,h (7) 

k 

Ck,a{l) = A k (l)£k,„ + ^2 B k:q (l)b q+k cl (8) 



where b q = b q (l = 0) and c k ^ a = c ktCr (l = 0). The initial I = values are 

(9) 



A k (0) = 1, B k<q (0) = 0, 



x q (o) = i, y M (o) = o. 

Unknown parameters Y qk (l), B kq {l) can be found directly from the flow equations (^,^) 

dB ktq (l) _ 



dl 

dY q , k (l) 
dl 



= -a q , k (l)v q , k {l)X q+k {l)A q {l) (10) 
= a k>q - k (l)v k>q - k (l)A k (l)A q - k (l). (11) 



There is no unambiguous way for determining the other missing coefficients A k (l), X q (l) because with the Ansatz 
^ ne ^ ow e 1 ua ti° ns still do not close. In what follows bellow we determine these parameters using the constraint 

W0,^0)] = <W, (is) 

{cfc, CT (0)4>'(0} =&k,k'&a,o'- (13) 

These statistical commutation/anticommutation relations yield 

A k {lf + J2 BkM ihl - h +q ) = 1, (14) 

x q (if + ]r y,Ai) 2 (i - h,i - u-k A ) = i. (15) 

k 

Here f kt a, b q denote the Fermi and Bose distribution functions. They appear after introducing the normal ordered 
forms for operators left of the commutation (|lj) and anticommutation (|lj) when substituting (0j|). 

We solved numerically a set of the coupled equations (|l^,|ll],[l4 15) simultaneously with the additional flow equations 
for e1(l),E k (l) and v ktq (l) given in the Ref. ||. Initially, for I = 0, we assumed: a) the one dimensional tight binding 
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dispersion eZ(l = 0) = — 2t cosfc, b) localized bosons E q (l = 0) = As and c) local exchange interaction Vk,q(l = 0) = v. 
For the sake of comparison with the previous studies of this model we used D = At = 1, A# = —0.6, v = 0.1 and 
total concentration of carries fixed to be 1. 
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FIG. 1. Momentum dependence of the fermion coefficient Afe(oo) for three temperatures T = 0.01, 0.02 and 0.03 (in units of 
the initial fermions bandwidth). Note a strong renormalization of the coefficient Ak around the characteristic momentum k* . 

In figure 1 we plot the momentum dependence of coefficient for I = oo (when fermion subsystem is decoupled 
from the boson one). Considerable variation of the parameter A^ from its initial value (being 1) can be observed 
for such momenta which are located near k* ~ 0.89. As discussed in the Ref. || this situation corresponds to the 
resonant scattering processes when Ek* + £-k* — It is worth mentioning that for ntot = 1 the Fermi momentum 
approaches this k* from bellow when temperature decreases. We could think of the coefficient Ak as a measure of 
the spectral weight for fermions. The missing part |1 — Ak\ is transfered towards some composite objects (mixtures 
of fermions and bosons or fermion pairs). Taking into account the fact that for a decreasing temperature kp shifts 
closer and closer to k* , we shall thus observe a diminishing spectral weight of fermions. This is in agreement with 
interpretation of the ARPES experimental results Q . 
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FIG. 2. Momentum dependence of the boson coefficients X g (oo) and Y^fc(oo) tor the same temperatures as in figure 1. 
The most efficient renormalization of boson coefficients occurs near q ~ 2k* . 

For a completeness we illustrate also the resulting {I = oo) values of the boson coefficients given in (0). Parameter 
A 9 shows only a small deviation from its initial value 1. This time such a change is most effective for boson momenta 
close to 2k*. It is exactly the same region of the Brillouin zone where the effective boson dispersion shows a kink (see 
figure 3). 

The missing part of the boson spectral weight is transferred to the fermion pairs, as given by the Ansatz (]?]). Bottom 
figure 2 shows that Y q ,k (which roughly measures the fermion pairs spectral weight) changes a lot for q ~ 2k*. 
Again, if we recall that for a decreasing temperature kp — > k* we can conclude that there is an increasing probability 
of finding more and more fermion pairs, even though our system is not in the superconducting phase. This is one of 
possible ways for observing the precursor effect. 
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FIG. 3. Effective boson dispersion Ek(l = oo) as measured from 2p level for the same set of temperatures as in figure 1. 
Recall that initially (for I — bosons were localized. 

In conclusion we studied evolution of fermion and boson operators in the BF model. During a continuous transfor- 
mation which eliminated the exchange interaction we observe transformation of the initial fermion and boson particles 
into the composite objects. In particular we find that fermion (boson) spectral weights at kp {2k f) are reduced when 
temperature decreases. This effect is in an agreement with the recent interpretation of the ARPES measurements for 
underdoped Bi2212 cuprates. To clarify a problem of the quasiparticles life time we need a more detailed study of the 
correlation functions. Such an analysis is currently under investigation and the results will be published elsewhere. 
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